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Bipolar ejection by the symbiotic binary system Z And 
during its 2006 outburst* 



o 
o 

o 

Q 

OO 



> 

O 
in 

(N 
O 

6 

c3 



N. A. Tomov 1 !, M. T. Tomova 1 and D. V. Bisikalo 2 

1 Institute of Astronomy, National Astronomical Observatory Rozhen, POBox 136, 4100 Smolyan, Bulgaria 
2 Institute of Astronomy of the Russian Academy of Science, 48 Pyatnitskaya Str., 119017 Moscow, Russia 



Accepted 2006 December 13. Received 2006 December 12; in original form 2006 November 24 



ABSTRACT 

High resolution data in the region of the line Ha have been obtained at the time of the 
light maximum and after it of the 2006 optical outburst of the symbiotic binary Z And. 
A blue-shifted absorption component indicating outflow velocity of about 1400 kms -1 
as well as additional emission components with similar velocities, situated on the two 
sides of the main peak of the line were observed during that time. It is suggested that 
all of them are spectral signature of bipolar outflow, observed for the first time in the 
optical spectrum of this binary. The emission measure and the mass of the nebular 
part of the streams were estimated approximately, which reached values of up to about 
6 xlO 57 cm" 3 and 8 xlO" 8 M . 

Key words: binaries: symbiotic - stars: mass-loss - stars: winds, outflows - stars: 
individual: Z And. 



1 INTRODUCTION 

Symbiotic stars are interpreted as interacting binaries con- 
sisting of a cool visual primary and a compact secondary 
component accreting matter from the atmosphere of its com- 
panion. As a result of accretion the compact object under- 
goes multiple outbursts accompanied with intensive mass 
ejection in the form of optically thick shells, stellar wind, 
single discreet ejections (blobs) as well as collimated bipo- 
lar jets. The system Z And is considered as a prototype of 
the classical symbiotic stars. I t consists of a normal c ool 
giant of spectral type M 4.5 (|Murset fc Schmidl [l999h . a 
hot compact compon ent with a temperature of 1.5x10 s K 
|Sokoloski et al . 2006) and an extended circumstellar nebula 
formed by the winds of the stellar components and partly 
photoionized by the compact object. 

Z And has underwent several active phases (after 
1915, 1939, 1960, 1984 and 2000) characterized by re- 
peated optical brightenings, partly realized as a result 
of a strong redistribution of the energy of the compact 
object. The energ y redistribution was caused by inten- 
sive loss of mass (ISwings fc Struvelll94ll; lBovarchuklll967l; 



Fernandez-Castro ct al 
2003a; ISokoloski et al 



.19951 ; iTomov. Taranova fc Tomoval 
2006J ; ISkopal et al.1 120061 ). The last 



active phase of Z And began in the end of August 2000 
and continues up to that time. The system underwent four 
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consecutive optical brightenings, whose maxima were in De- 
cember 2000 (V ~ 8.8 mag), November 2002 (V ~ 9.8 mag), 
Septe mber 2004 (V ~ 9 1 mag ) and July 2006 {V ~ 8.6 
mag) (jSkopal fc Pribullal [2006] . our UBV data). Spectral 
indications for intensive loss of mass were obtained dur- 
ing three of them, whose maxima were in 2000, 2002 and 
2006. The UV line Pv 1117 had a variable P Cyg pro- 
file containing one or two absorption components indicat- 
ing velocities of 0-300 kms" 1 in November and Decem- 
ber 2000. In all cases, however, the absorption reached 
its maximal depth in the velocity interval 0-50 kms 
(ISokoloski et al. I l2006h ■ During the same time the profiles 
of the lines He I were of type P Cyg too and the ab- 
sorption component i ndicated moderate velocities of the 
flow of 50-100 km s" 1 (|Tomov. Tomova fc Zamanovll2003"bl ; 
ISkopal et afl l2006h During the period September - Novem- 
ber 2002 the line Hell 4686, from its side, had two compo- 
nent emission profile consisting of a central narrow compo- 
nent and a broad component i ndicating stellar wind with a 
veloc it y of 1100-120 kms " 1 l|Tomov. Tomova fc Taranova! 
l2005h . ISkopal et al.l (|2006h concluded about intensive loss 
of mass by the compact object, realized at that time too, 
analyzing the wings of the line Ha. 

Our data show that the lines of H I and He I had P Cyg 
profile with absorption component indicating moderate ve- 
locity during the period July - September 2006. In some 
cases the components were more than one. In July, however, 
these lines had additional absorption component with a high 
velocity of about 1400 kms -1 as well. Moreover, during the 
period July - September the Ha line had additional emission 
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Table 1. Journal of observations. 



Date 


JD- 


Orb. 


Exp. Time 




2 453 000 


phase 


(min) 


Jul 8 


924.544 


0.837 


15+5 


Jul 9 


926.431 


0.840 


10+20 


Jul 14 


931.464 


0.846 


10+20 


Jul 19 


935.575 


0.852 


10 


Aug 8 


956.392 


0.879 


5+20+10 


Aug 12 


960.455 


0.885 


10+20+20+5 


Sept 7 


985.594 


0.918 


10+20 


Sept 8 


986.510 


0.919 


20 



components, placed on the two sides of its central peak and 
having velocities of 1200-1500 kms - . The appearance of 
the H/3 line was similar. The additional e mission line compo- 
nents of Z And were already repo rted bv lSkopal fc Pribullal 
(|2006h and ISkopal fc Wold (|200fJ ). which consider them as 
indications of bipolar jets. We suppose that all components 
indicating high velocities appear in bipolar outflow from the 
outbursting hot object. The 2006 outburst of Z And is a 
first one, when optical bipolar outflow from this system was 
observed. Our note is devoted to this phenomenon, where 
preliminary results of one study of the Ha line are treated. 



2 OBSERVATIONS AND REDUCTION 

The region of the line Ha was observed on eight nights in 
July, August and September 2006 with the Photometries 
CCD camera mounted on the Coude spectrograph of the 
2m RCC telescope of the National Astronomical Observa- 
tory Rozhen. During the last of these nights besides the 
region of the line Ha, that of H/3 was also observed. The 
spectral resolution was 0.2 A px on all occasions. Some 
of the exposures were comparatively long to obtain better 
ratio signal to noise of the continuum although the central 
emission component of the Ha line was saturated in some 
cases (Table 1). Ever when we made more than one exposure 
per night, the spectra were added with the aim to improve 
the signal to noise ratio. The IRAF package was used for 
data reduction. The additional emission components were 
fitted with a Gaussian to measure their radial velocity and 
equivalent width. It was made only for those spectra where 
they are intensive enough (Table 2). The inner uncertainty 
of the equivalent width in these cases is not more than 
50 per cent. The Ha flux was obtained using the equiva- 
lent width and the R c flux supposing that it is practically 
equal to the continuum flux at the positi on of this line. For 
this a im we used the R c magnitu des of ISkopal fc Pribullal 
(|2006h and ISkopal fc WoU i|2006t ) taken very close to the 
time of our observations. The uncertainty of the continuum 
flux is not more than 10 per cent and that of the line flux 
- not more than 50 per cent. The fluxes were corrected 
for an interstellar exti nction of E{B — V) = 0.30 accord- 
ing to the approach of ICardelli. Clayton fc Mathisl (|l989f ). 
It was used the ephemeris Min(vis) = JD 2 442 666 d + 
758*? 8 x E, where the orbital period is based on both pho- 
tometric and spectral data (iFormiggini fc Leibowitj Il994l; 
iMikolaiewska fc Kenvon1ll996l ; iFekel et al.ll2000T l. 
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Figure 1. Time evolution of the Ho line. 

3 ANALYSIS AND DISCUSSION 

In our spectra the line Ha had a strong central component 
located around the reference wavelength, broad component 
and additional absorption and emission features on both 
sides of the central component (Fig. [TJ. The central com- 
ponent was single-peaked having shoulder(s) on its short- 
wavelengths side, which was not visible only on the spectra 
taken in September. A weak peak component on the short- 
wavelengths side of the central component is seen on the 
spectrum of August 8 and the dip feature between them 
indicates a moderate velocity of about 100 kms -1 like the 
P Cyg absorption of the He lines. That is why this dip fea- 
ture can be determined by an outflow too, but it does not 
dispose below the level of the continuum probably because 
of filling by the emission. 

The broad component was not located around the refer- 
ence wavelength during all time of observations and in July 
its red wing was appreciably more intensive than the blue 
one. In our treatment, however, we will concentrate on the 
other additional components of the Ha line, which appeared 
during the period of our observations. 

The first spectra in July (Fig. [1} show a pronounced 
absorption with a velocity of —1400 kms -1 on the short- 
wavelengths side of the central Ha component and only a 
weak emission feature with irregular form and a velocity of 
about 1500 kms" 1 on its long- wavelengths side. We will con- 
sider the absorption first supposing that it is due to mass 
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flow from the compact object in the system. In this case 
the outflowing material generates an observed photosphere 
(pseudophotosphere, shell) whose optical depth in the con- 
tinuum is close to unity. If we assume that the mass flow 
has spherical symmetry and a constant velocity, this photo- 
sphere will be a spherical layer with an internal boundary R 
equal to the radius of the star and an external boundary of 
infinity. In this case the density is a function of the distance 
to the center and is expressed via the continuity equation. 
The column density N is expressed as 



2.0 



TV 



n(r) dr 



M 



where = 1.4 is a paramete r determining the mean mo lecu- 
lar weight fimii in the wind IjNussbaumer fc Vogellll987l ) and 
the other quantities have their commonly accepted meaning. 
The optical depth of the layer is r = kN , where k is the con- 
tinuum absorption coefficient (opacity). The mass- loss rate 



M = 



4tc^ijib_vRt 



Let us calculate the mass-loss rate. The opacity k per one 
atom is determined from the opacity krt per unit mass by 
means of the relation k — kR/im,H. We used a value of 0.5 
cm 2 g _1 of the Rosseland mea n opacity in the visu al region 
from Table 4 from the work of ISeaton et al.l l| 19941 ). In this 
case k — 1.2 x 10~ 24 cm 2 is obtained. However, we have 
no estimate of the size of the observed photosphere of the 
hot component of Z And at the time of our spectral obser- 
vations. A radius of the hot co mponent of 2.4(d/1.12 kpc) 
Rq was obtained in the work of iTomov et al.l (2003a) from 
analysis of the continuum energy distribution of the system 
at the time of the light maximum of its 2000 outburst. As a 
result of the analys is of the spectral energy distribution too 
ISkopal et al.l (|2006h obtaine d a radius of 3 . 1(d/1.5 kpc) Rq. 
The correction of the data of lTomov et all 11200331 ) for a dis- 
tance to the system of 1.5 kpc leads to value o f 3.2 R^ and 
is in v ery good agreement with the result of ISkopal et alj 
(2006). For the need of our approximate calculation we will 
use the value of 3.5 Rq, as according to our UBV data the 
2006 outburst is more luminous than the 2000 one, when 
the increase of the l ight was due mostly to the expansion of 
the hot component (|Tomov et al.ll20 03a|) . Having all of these 
quantities, for the mass-loss rate we obtain M=1.33xl0~ 5 
Mq yr- 1 . 

This value is too great compared with the ob- 
served mass-loss rate of the hot compact component 
of the symbiotic systems during their active phases, 
which is most frequently about 10~ 9 - 10~ 7 M^yr" 1 
dVogel fc Nussbaumei] Il994l; lAftamore fe Cassatellal1l997l ; 
iGuoliang. Yungelson fc H an 2006). The great value follows 
from our assumption for spherical symmetry of the outflow- 
ing material and the observed high velocity. It can be in 
reasonable limits if the material does not flow in all direc- 
tions, but only in some of them, related to areas with small 
angular sizes. The material can flow mainly in the polar ar- 
eas, if it is prevented in equatorial one for example by an 
accretion disk. In this case the velocity in the polar areas in- 
creases to provide the mass ejection there. This supposition 
can be related to the existence of two kinds of absorption 
components indicating very different expansion velocities. 
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Figure 2. Transition from the absorption feature to emission one 
indicating the front stream. The level of the local continuum is 
marked with a dashed line. 



We can assume that the dip feature of Ha is due to mass 
flow in the orbital plane with a moderate velocity. The flow- 
ing material, however, is prevented by the accretion disk and 
the direction of its motion changes. In this way it is ejected 
in a direction, perpendicular to the plane of the disk, like 
the gas flowing in that direction. All of the material gives 
rise to the high velocity components of the Ha line. Thus 
the blue shifted absorption indicates motion probably in the 
front part of the wind which is projected on the stellar disk. 
Then we can suppose that the weak emission feature with a 
velocity of 1500 kms -1 appears in the nonocculted part of 
the back component of the wind. 

As it is seen from the evolution of the spectrum in 
Figs. [T] and [2] the blue-shifted absorption component disap- 
pears and an emission rises. Thus two emission components 
on the two sides of the central peak form in July. The dis- 
appearance of the blue-shifted absorption and the rise of an 
emission is most probably due to decrease of the flow rate 
and/or increase of the number of the emitting atoms in that 
part of the wind which does not project on the stellar disk. 
If the absorption is related to the inner part of the wind and 
the emission - to its outer one, not projecting on the stellar 
disk, the evolution of the spectrum shows that the velocity 
in the front wind component (stream) decreases with the 
distance to the star. 

It is seen in Fig. [1] that the emission on the long- 
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Table 2. Radial velocities, equivalent widths, line fluxes and emission measures of the outflow components of the 
Ha line. 









Blue 








Red 




Date 


RV 


EW 


FxlO -12 


n 2 V x 10 57 


RV 


EW 


FxlO -12 


n 2 V x 10 57 




(kms" 1 ) 


(A) 


(erg cm -2 s — 1 ) 


(cm -3 ) 


(km s —1 ) 


(A) 


(erg cm -2 s _1 ) 


(cm -3 ) 


Jul 19 


-1196 


0.8 


1.754 


2.4 


1445 


0.9 


1.813 


2.5 


Aug 8 


-1087 


1.5 


2.912 


3.9 


1346 


1.8 


3.442 


4.7 


Aug 12 


-1260 


2.5 


4.828 


6.5 


1201 


2.0 


3.846 


5.2 


Sept 7 


-1245 


1.8 


3.109 


4.2 


1178 


2.0 


3.362 


4.6 


Sept 8 


-1262 


1.7 


2.881 


3.9 


1196 


2.1 


3.534 


4.8 



wavelengths side of the central peak increases compared to 
the continuum. It is seen also in the figure that the veloc- 
ity of the line decreases from about 1500 kms -1 to about 
1200 kms - . Taking into account the fact that the distance 
covered by the emitting gas increases with time it can be 
concluded that the velocity in the back stream decreases 
with the distance to the star too. 

The H/3 line had similar emission components, as it is 
seen on the unique frame taken in its region in September. 

To know the energy emitted in the Ha high velocity 
components as well as the quantity of the emitting gas we 
calculated their line flux and emission measure. The emis- 
sion measure was obtained supposing that the gas medium is 
optically thin and has a constant density. It is not possible 
to obtain the electron temperature and the electron den- 
sity in the region where the outflow components are emitted 
from observation since we have no indication about the ap- 
pearance in this region of certain lines giving information 
for these parameters. We suppose that the electron temper- 
ature is 20 000 K and the mean electron density is about 



10" 



cm in the nebular portion of the streams. 



Then 

.3 „-l 



we used recombination coefficient of 5.956x10 
for case B, corresponding to these temperature and den- 



sity (|Storev fe H ummer 1995) . We adopt ed a helium abun- 
dance of 0.1 ( Vogel fe Nussbaumei Il994 | ) and a distance to 
the sy stem of 1.5 kpc according to Mikolaiewska fe Kenvonl 
(1996). For the calculation of the emission measure we need 
to know the state of ionization of helium in the emitting 
region. The line Hen 4686 was very weak (as it is seen from 
our data) during the period of our observations which means 
that singly ionized helium is dominant in the circumbinary 
nebula at the times close to the epoch of the maximal light. 
That is why we assume the state of ionization to be He + in 
the region of the streams. 

The line flux and the emission measure are listed in Ta- 
ble 2. These data propose that the ejected matter increases 
till August 12 and after that decreases, which is probably 
due to decrease of the flow rate. We calculated also the 
mass of the nebular portion of the streams on the basis of 
the greatest emission measure observed on August 12. It 
turned out to be about 8xl0 -8 Mq for the front stream 
and 6xl0 -8 Mq for the back one. 



4 CONCLUSIONS 

We present results of high resolution spectral observations 
carried out in the region of the line Ha of the spectrum 
of the symbiotic binary Z And close to the light maximum 



of its 2006 optical brightening. The line Ha had a blue- 
shifted absorption component with a velocity of about 1400 
kms -1 at the time of maximal light in July. This line had 
also additional emission components placed at 1200-1500 
kms -1 on the two sides of its central peak in July, August 
and September. The line H/3 had similar components as far 
as it is seen on its unique frame. It is supposed that both 
the absorption and the emission components are signature 
of bipolar wind outflow from the hot compact object in this 
system. The different velocities observed at different spectra 
indicate probably the decrease of the flow velocity. 

The energy flux of the outflow components of the Ha 
line and their emission measure were calculated on the basis 
of -R c photometric data from the literature taken during the 
time of the spectral observations. These quantities increased 
till August 12 and after that decreased, which is probably 
due to diminution of the flow rate. The mass of the nebu- 
lar part of the streams on August 12 was obtained, which 
amounts to 6-8 xlO" 8 Mq. 
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